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Abstract. Random motions can occur in the intergalactic gas of galaxy clusters at all stages of their evolution. Depending 
on the poorly known value of the Reynolds number, these motions can or cannot become turbulent, but in any case they 
can generate random magnetic fields via dynamo action. We argue that magnetic fields inferred observationally for the intra¬ 
cluster medium require dynamo action, and then estimate parameters of random flows and magnetic fields at various stages 
of the cluster evolution. Polarization in cluster radio halos predicted by the model would be detectable with the SKA. 

Key words: Galaxies: clusters; general - galaxies: intergalactic medium - galaxies; magnetic fields - MHD - turbulence 
©0000 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 


1. Introduction 

We witness the epoch when galaxy clusters are being formed. 
Galaxies and smaller-size structures have generally achieved 
(quasi-)steady states in their evolution, whereas structures of 
larger scales keep evolving. There is abundant evidence of 
formation processes in galaxy clusters, including the merger 
of large masses comparable to the cluster mass. Relaxed, 
symmetric galaxy clusters are rare. Thus, most constituents 
of the clusters keep evolving. Our particular interest here is 
with motions of the gas and magnetic fields in the intergalac¬ 
tic space of evolving galaxy clusters. The epoch of major 
mergers is notable for widespread, intense random motions, 
followed by a period of decaying flows. In a steady state, ran¬ 
dom motions can be confined to the wakes produced by the 
cluster galaxies and smaller mass subcluster clumps that con¬ 
tinue falling into the cluster; as we argue here, these motions 
are likely to occupy a small fraction of the total volume. 

The high electrical conductivity of the intracluster gas is 
not sufficient to ensure the survival of a magnetic held cap¬ 
tured by the forming cluster. The reason is that the intracluster 
gas is very viscous, and any inhomogeneous magnetic held 
would drive gas motions which will rapidly decay, thereby 
converting magnetic energy into heat. Even for a small vis¬ 
cosity, the held would decay by driving decaying MHD tur¬ 
bulence. On the other hand, random gas motions of the intra- 
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cluster gas can be efficient generators of (random) magnetic 
fields via a process known as the fluctuation dynamo. 

The origin and properties of magnetic held in the intra¬ 
cluster plasma cannot be understood without knowing the na¬ 
ture and parameters of the plasma motions. The reason is that 
the magnetic Reynolds number in the intracluster gas is large, 
hence magnetic induction effects are strong for any plausible 
speed of gas motion. Our model of gas motions and mag¬ 
netic fields in galaxy clusters, presented and justified in de¬ 
tail by Subramanian et al. (2006), is in close agreement with 
the available data. It predicts that magnetic held in the intra¬ 
cluster gas is represented by magnetic sheets of a thickness 
20-30 kpc, with magnetic held strength about 2-4 fj,G within 
the sheets. These structures All about 20% of the volume of a 
single turbulent cell. Apart from this intermittent component, 
magnetic held also has a weaker, widely distributed part. 

We propose that, at late stages of cluster evolution, turbu¬ 
lence is confined to the wakes of galaxies and infalling mass 
clumps, and does not All the volume. However, the area cov¬ 
ering factor of the wakes can be close to unity, so most lines 
of sight pass through a turbulent region. This can help to rec¬ 
oncile the available evidence of turbulence in such clusters as 
Perseus with the existence of long filaments in the intracluster 
gas that seem to be inconsistent with pervasive turbulence. 

Since typical lines of sight through a cluster pass through 
a small number of turbulent cells, synchrotron emission from 
galaxy clusters can show a detectable degree of polarization 
at wavelengths 3-6 cm. 
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Fig. 1. Snapshots of magnetic field in a cross-section through the middle of the computational domain in a numerical sim¬ 
ulation of turbulence driven by an imposed random force and its dynamo action. The left-hand panel shows a statistically 
steady state at a time f/foi = 0.30 whereas the right-hand panel illustrates magnetic field structures in turbulence at a late 
stage of decay, f/foi ~ 60. Here tgi is the eddy turnover time before the start of the turbulence decay (given in Table^J. The 
dimensionless energy injection scale in these simulations is about 4 (with the domain size of 27r), so each frame contains a few 
turbulent cells. The strongest magnetic field within the frame is close to the equipartition value with respect to the turbulent 
energy. The magnitude of the field component perpendicular to the plane of the figure is shown color coded (in shades of 
grey) with black corresponding to field pointing into the figure plane, and lighter shades, to field pointing out of the plane. 
The field in the plane of the figure is shown with vectors whose length is proportional to the field strength. 


2. Three evolutionary stages 

2.1. The epoch of major mergers 

Theories of hierarchical structure formation suggest that clus¬ 
ters of galaxies have been assembled relatively recently. N- 
body simulations indicate that the clusters form at the inter¬ 
section of dark matter filaments in the large-scale structure, 
and result from both major mergers of objects of comparable 
mass (of order IO^^Mq) and the accretion of smaller clumps 
onto massive protoclusters. It is likely that intense random 
vortical flows, if not turbulence, are produced in the merger 
events. Their plausible properties are summarized in Table [2 
The structure of magnetic field at this stage is illustrated in 
the left-hand panel of Fig. [Q What is shown is the statisti¬ 
cally steady state of magnetic field produced by dynamo ac¬ 
tion in turbulent flow with the Reynolds number about 400 
and the magnetic Prandtl number equal to unity. Similar mag¬ 
netic structure plausibly occur in the turbulent wakes of sub¬ 
clusters and galaxies as well. 

It is not quite clear whether random flows driven during 
major merger events and at later stages of evolution will de¬ 
velop into turbulence. The nature of the flow depends on the 
value the Reynolds number which is difficult to estimate re¬ 
liably for the collisionless, magnetized plasma of the intra¬ 
cluster space where plasma instabilities can be responsible 
for anomalous viscosity and resistivity (Schekochihin et al. 
2005). The problem is further complicated by the possibility 
of dynamo action, since the magnetic field can affect both vis¬ 
cosity and magnetic diffusivity. This may lead to the growth 
of the magnetic diffusivity and reduction of viscosity as the 


magnetic field is being amplified by the dynamo, so that the 
magnetic Prandtl number tends to unity. 


2.2. Decaying turbulence 

Random flows produced by major mergers decay after the 
end of the merger event. Unlike a laminar flow that decays 
exponentially in time due to viscosity, turbulent kinetic en¬ 
ergy decays slower, as a power law (e.g., Landau & Lifshitz 
1975; Frisch 1995). The reason for this is that kinetic energy 
mainly decays at small scales, to where it is constantly sup¬ 
plied by the turbulent cascade. As a result, the energy decay 
rate depends nonlinearly on the energy itself, which makes 
the decay a power law in time. Our simulations confirm that 
the power-law decay occurs even for the Reynolds number as 
small as Re « 100. At this stage of evolution, the turbulent 
scale Iq grows with time, whereas turbulent energy density E 
reduces, together with the turbulent speed vq, typically as 

F; ~ iup oc lo oc for f/foi » 1 , 

where subscript ‘i’ refers to the start of the evolution, foi is a 
certain dynamical time scale, which can be identified with the 
initial turnover time of the energy-containing eddies, tgi = 
loi/voi, subscript ‘0’ refers to the energy-range (correlation) 
scale of the motion. The structure of magnetic field in the 
decaying flow is shown in the right-hand panel of Fig.Q and 
parameters of the flow and magnetic field are shown in the 
second line of Tabled] 
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Table 1. Summary of turbulence and magnetic field parameters at various stages of cluster evolution: duration of the stage (the 
last two stages represent steady states), the r.m.s. velocity vo and scale Iq of turbulence and eddy turnover time to = Iq/vq (for 
the decaying turbulence, values for the middle of the decay stage are given, 2 Gyr after its start), the equipartition magnetic 
field Beq = ( 47 rpuQ)^/^ with p the gas density (i.e., maximum field strength within a turbulent cell), thickness of magnetic 
filaments and sheets Ib for the statistically steady state of the dynamo, the r.m.s. magnetic field within a turbulent cell Srms, 
and finally the standard deviation of the Faraday rotation measure ctrm (calculated for the volume filling turbulence along 
path length of 750 kpc through the central parts of a cluster in the first two lines, and assuming one transverse wake along the 
line of sight in the last two lines). A subcluster mass of 3 x has been assumed. 


Evolution stage 

Duration 

(Gyr) 

Vo 

(kms“^) 

^0 

(kpc) 

^0 

(Gyr) 

Beq 

(pG) 

Ib 

(kpc) 

-^rms 

(f^G) 

CRM 

(radm~^) 

Major mergers 

4 

300 

150 

0.5 

4 

25 

1.8 

200 

Decaying turbulence 

5 

130 

260 

2.0 

2 

44 

0.8 

120 

Subcluster wakes 


260 

200 

0.8 

4 

34 

1.6 

110 

Galactic wakes 


300 

8 

0.03 

4 

1.4 

1.6 

5 


2.3. Tfirbulent wakes of subclusters and galaxies 


At the final stage of the evolution, when the cluster enters a 
steady state, turbulence is maintained only in the wakes of 
galaxies and smaller mass clumps that continue to accrete 
onto the cluster. The wakes become weaker as the gas within 
the clumps or galaxies is stripped by the ram pressure of intra¬ 
cluster gas. The radius of a wake at its head is close to the 
radius within which gas of the mass clump or galaxy remains 
intact. We estimate the stripping radius as Rq ~ 100 kpc 
for clumps of a mass 1 O^^M 0 (which fall into a cluster ev¬ 
ery 3 Gyr) and Rq = 3-5 kpc for massive elliptical galax¬ 
ies. If the flow within the wake becomes turbulent (so that 
it can be described in terms of Prandtl’s theory of turbulent 
wakes), the wake length X is controlled by the magnitude of 
the Reynolds number via 

X/Rq ~ (Rei/Rccr)^ , 


where Rccr ~ 400 (Tomboulides & Orszag 2000) is the 
marginal Reynolds number with respect to the onset of tur¬ 
bulence. This value of Rccr was obtained for a flow around a 
solid sphere; Rccr for gas spheres is not known. The strong 
dependence of the wake parameters on the Reynolds number 
makes the estimates somewhat uncertain. On the other hand, 
it implies that galactic wakes can be very sensitive to the de¬ 
tailed parameters of the galactic motion and intergalactic gas, 
so that clusters with very similar parameters can have vastly 
different wake structures. 

The area covering and volume filling factors, fs and /y, 
respectively, of W = 5 wakes, produced by 10^^Mq subclus¬ 
ters, within the virial radius r « 3 Mpc, are estimated as 

N ^ Ro V ^RecA~‘^ f X 


/s-0.15- 


5 \ 100 kpc 


400 


1 kpc 


N 

fv ^ 0.02 - 
5 


Rq 


100 kpc 


RCcr 

loo 


-5 


-5 


Ikpc 


where W « 5 is consistent with models of hierarchical 
structure formation, and A is an effective mean free path in 
the intracluster gas (its introduction is an attempt to allow 
for our insufficient understanding of viscosity mechanisms). 
The covering and filling factors strongly depend on A and 


RCcr- Furthermore, both fs and fy depend on high powers 
of another poorly known parameter, the stripping radius Rq. 
Hence, as noted above, properties of the subcluster wakes can 
be rather different in apparently similar clusters. In addition, 
numerical simulations of turbulent wakes should be treated 
with caution as otherwise reasonable approximations, numer¬ 
ical resolution, and numerical viscosities can strongly affect 
the results. On the other hand, it is plausible that fs = 0(1) 
but fv 1 , so that a typical line of sight through the cluster 
intersects at least one turbulent region (where our estimate of 
the r.m.s. turbulent speed is 200-300 km s“^) despite the fact 
that turbulence occurs only in a small fraction of the clus¬ 
ter volume. Thus the presence of long Ha filaments observed 
by Fabian et al. (2003, 2006) in the core of the Perseus clus¬ 
ter may not be inconsistent with various evidence for random 
motions in this cluster core (Churazov et al. 2004; Rebusco et 
al. 2006). A possible signature of such spatially intermittent 
turbulence could be a specific shape of spectral lines, with a 
narrow core, produced in quiescent regions, accompanied by 
nonthermally broadened wings. 

The area covering factor of galactic wakes within the gas 
core radius, 180 kpc, is unity if 

X/Ro ~ 30-15 , a: ~ 100-70 kpc , (1) 

and the volume filling factor of such wakes is fv — 0.07. 
The length of galactic wakes required to cover the projected 
cluster area, given by Eq. Q, does not seem to be unre¬ 
alistic. For example, Sakelliou et al. (2005) have observed 
a wake behind a massive elliptical galaxy (mass of order 
2 X lO^^M©) moving through the intracluster gas at a speed 
about Vc — 1000 km s“^. The length of the detectable wake 
is about X ~ 130 kpc (assuming that it lies in the sky plane), 
and its mean radius is 40 kpc (obtained from the quoted vol¬ 
ume of about 2 X 10® kpc®). The projected area of the wake 
is about 10 "^ kpc^, as compared to 10 ® kpc^ for the wake pa¬ 
rameters derived above. This wake has been detected only 
because it is exceptionally strong, and it is not implausible 
that weaker but more numerous galactic wakes can cover the 
projected area of the central parts of galaxy clusters. 

We conclude that subcluster wakes are likely to be tur¬ 
bulent, but galactic wakes can be laminar if the viscosity of 
the intracluster gas is as large as Spitzer’s value. Given the 
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uncertainty of the physical nature (and hence, estimates) of 
the viscosity of the magnetized intracluster plasma, we sug¬ 
gest that turbulent galactic wakes remain a viable possibility. 
Both types of wake have low volume filling factor but can 
have an area covering factor of order unity. Parameters of tur¬ 
bulence and magnetic fields produced within the wakes are 
given in the last two lines of Tabled 

3. Intracluster magnetic fields and their 
observational signatures 

A random flow of electrically conducting fluid (either turbu¬ 
lent or not) can generate and maintain magnetic field if the 
magnetic Reynolds number i?ni = volo/rj exceeds a cer¬ 
tain critical value, > iJm.cr where iJm.cr depends on the 
statistical properties of the flow and usually remains within 
the range i?ni,cr = 30-100 (here p is the magnetic diffusiv- 
ity). This condition is undoubtedly satisfied in the intracluster 
plasma, where ^ 1 (unlike the kinematic Reynolds num¬ 
ber Re which may be quite modest), and so the dynamo ac¬ 
tion in galaxy clusters is more than plausible. This type of dy¬ 
namo, known as the fluctuation dynamo because it produces 
fluctuating magnetic fields with zero mean value, does not re¬ 
quire any overall rotation, density stratification, a-effect, etc. 
Reviews of fluctuation dynamos can be found in Zeldovich et 
al. (1990) and Brandenburg & Subramanian (2005). 

At an early stage of the dynamo action, when magnetic 
fields are still too weak to affect the flow, magnetic energy 
density grows exponentially. The e-folding time of the r.m.s. 
magnetic held generated by motions of a scale I and speed v 
is of the order of t ~ Z /v. In Kolmogorov turbulence (and in 
any random flow with a sufficiently steep energy spectrum), 
the e-folding time is shorter at smaller scales, so that mag¬ 
netic energy spectrum has maximum at small scales, close 
to the magnetic diffusion scale = IqR^ with d = 1/2 
(Kazantsev 1967; Zeldovich et al. 1990; Schekochihin et al. 
2004) or d = 3/4 for the marginal mode in Kolmogorov tur¬ 
bulence (Brandenburg & Subramanian 2005 and references 
therein). Thus, magnetic held is concentrated into magnetic 
sheets (and filaments) of thickness Irj and length of order Iq. 
The held strength within the magnetic structures is, at satura¬ 
tion, close to equipartition with kinetic energy, B ~ i?eq = 
(dTrpvg)^/^, where p is the gas density. An application of 
this theory to galaxy clusters is discussed by Ruzmaikin et 
al. (1989), although these authors presumed that turbulence 
produced in galactic wakes can be volume-filling. 

The nonlinear stage of the fluctuation dynamo, where the 
growth of magnetic energy saturates because of the action of 
the Lorentz force on the flow, is more controversial (Branden¬ 
burg & Subramanian 2005). We consider plausible a model 
of the nonlinear dynamo suggested by Subramanian (1999); 
this model is consistent with numerical simulations of the 
fluctuation dynamo, especially for systems with unit mag¬ 
netic Prandtl number (Subramanian et al. 2006 and references 
therein). In this model, the thickness of magnetic structures 
in the saturated dynamo is of the order of the magnetic diffu¬ 


sion scale, but now based on the critical value of the magnetic 
Reynolds number, 

Ib — loRm]Ir ■ 

This scale is given in Table [^for various stages of evolution, 
together with the magnetic held within magnetic structures, 
Beq and the resulting r.m.s. magnetic held within a turbu¬ 
lence cell, i?rms- Table^also contains estimates of the r.m.s. 
value of the Faraday rotation measure, ctrm produced in the 
dynamo-generated magnetic held; the estimates agree with 
observations fairly well. Our numerical simulations of the 
fluctuation dynamo, where we directly measure urm, also 
confirm our semi-analytical estimates. 

Our model implies that the correlation scale of random 
motions in the intracluster gas, Iq, is larger than that assumed 
earlier. With Iq ~ 150 kpc (Table [0, only 5 turbulent cells 
occur along a path length of L = 750 kpc. The resulting de¬ 
gree of polarization of radio emission from clusters with syn¬ 
chrotron halos can be estimated as p ~ ~ 0.2, 

where po Ki 0.7, n ~ L/Iq is the number of magnetic 
structures along the line of sight (assuming that one mag¬ 
netic sheet with well-ordered magnetic held occurs in each 
turbulent cell and that the linear resolution is better than 
Iq), and a factor 1/2 allows, in a very approximate man¬ 
ner, for the volume-filling magnetic held outside the mag¬ 
netic sheet which only produces unpolarized emission. De¬ 
polarization by Faraday dispersion and beam depolarization 
can reduce the degree of polarization to a fraction of percent 
at long wavelengths. However, polarization observations at 
wavelengths 3-6 cm (where Faraday depolarization is suffi¬ 
ciently weak) can reveal magnetic structures produced by the 
dynamo action if the angular resolution is high enough. Such 
observations can become feasible for many galaxy clusters 
with the advent of the SKA (Feretti & Johnston-Hollitt 2004). 
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